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Abstract. We investigate a two-dimensional statistical model of A'^ charged particles 
interacting via logarithmic repulsion in the presence of an oppositely charged compact 
region K whose charge density is determined by its equilibrium potential at an inverse 
temperature corresponding to /3 = 2. When the charge on the region, s, is greater 
than A'^, the particles accumulate in a neighborhood of the boundary of K, and form a 
determinantal point process on the complex plane. We investigate the scaling limit, as 
Af — > oo, of the associated kernel in the neighborhood of a point on the boundary under 
the assumption that the boundary is sufficiently smooth. We find that the limiting kernel 
depends on the limiting value of N/s, and prove universality for these kernels. That is, 
we show that, the scaled kernel in a neighborhood of a point Q £ dK can be succinctly 
expressed in terms of the scaled kernel for the closed unit disk, and the exterior conformal 
map which carries the complement of K to the complement of the closed unit disk. When 
N/s — > we recover the universal kernel discovered by Lubinsky in [13]. 



1. Introduction 

1.1. Potential Theoretic Weights. Let iiT C C be a compact subset whose boundary 
T = dK is a Jordan curve. We will assume that T is sufficiently nice in a way that will be 
made precise in the sequel. For such K, there exists a unique measure ujk, the equilibrium 
measure on that minimizes the energy functional /[cr] := — JJloglz — u\da{u)da{z) 
among all positive probabihty measures a supported on K [17]. We define Pk : C — (0, oo), 
by 



(1) Pk(2:) := cxp + j \og\z - u\dujK{u) 

This function is simply the rescaled exponentiated equilibrium potential of K . 

It can be verified that loj^ is supported on T, is identically one on K and, as z — )> c», 
Pk{z)/\z\ — >■ 7^^, where '■= exp{— /[w/^-]} is the logarithmic capacity of K. 

In this paper we will be interested in random vectors whose joint density is given by 

1 r ^ 1 

(2) f^A^CA) — I []7«(A„)| [] |A„-A„p; A e C^, 

^ 71—1 m<n 

where 

(3) w{X)^PKi\r^\ Zn:^[ [\\w{Xn)\\\\\n-Xm?dA''{X), 

and s > N (that is s is sufficiently large to guarantee that Zn is finite). Here and throughout, 
A and are Lebesgue measure on C and respectively. 
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We will often refer to the components of such random vectors as eigenvalues, since the 
joint density (2) can be thought of as a modification of the joint eigenvalue density of the 
ensemble of matrices with i.i.d. complex Gaussian entries. The eigenvalues of this latter 
matrix ensemble, originally introduced by Ginibre [10], have joint density given by (2) where 

In Section 3 we will give (i) a matrix model whose joint density of eigenvalues is given by 
(2) with weight given as in (3) for K equal to the closed unit disk, as well as (ii) models for 
more general K where the components of A represent the positions of electrostatic particles 
confined to the plane and in the presence of a field determined by K, and (iii) an ensemble 
of random polynomials chosen with respect to a height function determined by K whose 
roots are distributed as in (2). 

Our primary goal is to demonstrate that, in the double scaling limit as s and N approach 
infinity, the local statistics of the eigenvalues near a point on the boundary of K depend 
only on the limiting ratio of s and N, but are essentially independent of the specifics of K. 
This will follow from the asymptotic behavior of the reproducing kernel of L'^(w), which 
in turn follows from the asymptotics of the leading coefficient of the related orthonormal 
polynomials. When s = oo, w is simply the characteristic function of K and our results 
collapse to those given by Lubinsky [13] for the universality of reproducing kernels formed 
with respect to Bergman polynomials for K. 

1.2. Eigenvalue Statistics. We briefiy review some basic concepts for solvable ensembles 
of random matrices and how they relate to eigenvalue statistics. In this section we will 
assume that the joint density of eigenvalues is given by (2) where, for the purposes of this 
section, w : C — >■ [0, oo) is any non-negative function such that < Zjv < oo. 

We will suppose that S = {^1,^2, ■ ■ ■ ,^n} C C is a random set corresponding to the 
eigenvalues of a random matrix from our ensemble. (Or, what amounts to the same thing, 
S is the set corresponding to a random vector sampled from the density ^n)- Given a set 
E Q C we may construct a random variable X given by the cardinality of S n Given 
disjoint subsets Ei, E2, . ■ ■ , En we will let Xi, X2, ■ ■ ■ , Xn be the corresponding random 
variables. The n-th correlation function of our ensemble is defined to be i?„ : C" — >■ [0,oo), 
where 

E[XiX2---Xn]:^ I I ■■■ I i?„(A)dA"(A). 

It is straightforward to see that i?jv ^ NlflN. A less obvious exercise is to show that for 
< n < iV, 

(4) i?„ (A) = — ^ / i?^ (A V x) dA^-" (x) , 

(N - ny. JcN-,^ 

where A V x = (Ai, • • • , A„, xi, . . . , xn-u)- Many probabilities of interest can be expressed 
in terms of correlation functions. One particularly important example is the gap probability 
that there are no eigenvalues in E, 

(5) Prob{X = 0} = ^ I i?„(A) dA"(A). 

Equations (4) and (5) are valid for a wide variety of symmetric measures on C^. However, 
the presence of the square of the Vandermonde determinant which appears in (2) leads to 
additional structure which may be exploited. 
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Suppose ttq, TTi, . . 
That is, 



, , 7rjv_i are the orthonormal polynomials with respect to the weight w. 



7i'„(2)7r™(z)w(z) dA = S„ 



where, as usual, 6nm is 1 or depending on whether or not ri = m. The kernel of the 
ensemble is defined by 

N-l 

Kn{z,u) ^/w{z)w{u) ^ 7r„(z)7r„(u). 



n=0 



(Following Lubinsky's notation, we will reserve the symbol K]\[ for the unweighted analog 
of this kernel). In a celebrated result, Mehta and Gaudin [14] were able to express the 
correlation functions of ensembles with eigenvalue density (2) in terms of determinants of 
matrices formed from this kernel, 



(6) i?„(A) =det [i^Ar(A„Afc) 

(See also [22] for a more modern derivation). 



1.3. Universality. When N is large we expect that, with high probability, the eigenvalues 
will accumulate in a neighborhood of dK. Slightly more precisely, if C G 9K, then the 
number of eigenvalues in a disk of (small) radius e about ^ is proportional to N; the constant 
of proportionality is given by the integral of the equilibrium measure over the arc of dK 
contained in the disk. The exact details of this phenomenon will be explored in a subsequent 
paper, for now we use this only as intuition to guess the proper scale on which we expect 
Kpf to converge. 

From (5) and (6), the probability that there are no eigenvalues in a disk of radius e 
centered at ( is given by 

N 



(7) 



n=0 



det 



e^J^Ar,.(C + eAj-,C + eAfc) 



dA"{\} 



where D is the disk of radius 1 centered at the origin. Here we have made explicit that the 
kernel is dependent on s as well as N. 

Under the assumption that there are 0{N) eigenvalues in a neighborhood of then 
we should scale e like 1/iV in order for (7) to approach a non-trivial limit. That is, the 
limiting gap probability of there being no eigenvalues in a shrinking neighborhood with 
radius e = 1/A^ is given by 

N 



(8) 



lim > 



(-1)" 



det 



1 

7V2 



K 



N.! 



N 



dA"(A). 



Since s > N, this limit also depends on how s scales with N and we will assume that N/s 
converges to some £ e [0, 1]. 

If it can be shown that there is some limiting kernel Hq £ so that 

1 ~ / z u\ ~ 



7V2 

uniformly on compact subsets of 

(-1)" 



then (8) converges to 



oo 



n=0 



det 



].k=l 



dA"(A). 
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(See for instance [1, §3.4]). Our primary result here is that H,^^i exists, and is dependent on 
C and K in only the most trivial manner. More specifically, we will express -ff^/ in terms 
of the limiting kernel for the ensemble formed from the closed unit disk and the value of a 
conformal map from C \ if to C \ D evaluated at C. This is what is called universality for 
potential theoretic ensembles. 

We will also demonstrate that ii^,f is a convex combination of i?(^.o (Lubinsky's limiting 
kernel) and 

1.4. Potential Theoretic Orthogonal Polynomials. We denote the orthonormal poly- 
nomials for the weight P^^^ , s > 1, by {7r„,s}|jLo^^ • That is, these are polynomials with 
positive leading coefficients that satisfy 

(9) / 7r„_s(z)7r„_5(z)F^^*(z) = (5„,m. 

Jc 

The reproducing kernel for this system of polynomials is given by 

(10) KnA^'^-^) — 7r„,,(z)7r„,,(w), < [s - Ij , 
with the weighted kernel given by 

(11) Kn^z.u) Pj^'{z)Pj,'{u)KnAz,u). 

Our derivation of Hq j will follow from the asymptotics of -ffjv.s, which in turn will follow from 
the asymptotics of the orthogonal polynomials. These latter asymptotics are of independent 
interest, and they provide the other primary results of the paper. 



2. Statement of Results 

In what follows, we assume that T :— dK is a rectifiable Jordan curve which is either 
analytic or of class (7^+^'", where p is a nonnegative integer and a £ (0, 1). That is, the 
arclength function of T is p times continuously differentiable as a periodic function on the 
real line and its p-th derivative is a-H61der continuous. Denote by $ the conformal map of 
O := C\K onto O := C \ D such that $(oo) = oo and $'(cx)) > 0. In the case where T is 
an analytic Jordan curve we denote by p{T) < 1 a number such that has a univalent 
extension into \w\ > p{T). Moreover, we put Op :— <i>^^({|fi;| > p}) for each p > p{T). 

It is known that |$| is identically equal to Pk on O and therefore $'(cx)) = 1k^- Hence, 
orthogonality relations (9) can be rewritten as 

(12) / ■7Tn,s{z)Trm,s{z)dA + / 7r„,5 (z)7r„,5 (z) | $(z) | "^"dA = 5„™ , 

JD JO 

n,m < [s — 2J , where D is the interior domain of T. Since |$| > 1 in O, we can formally 
set 7r„.oo to be polynomials satisfying 

^n,oo (^)'^m,oo (^) ^nm- 

In a sense, potential theoretic polynomials 7r„_s can be considered as perturbations of 7r„.oo- 
The latter were initially studied by Carleman [2] who derived their exterior asymptotics 
(asymptotics in O) for the case T being an analytic Jordan curve. The results in [2] were 
subsequently extended by Suetin [21] to include C^+i'" Jordan curves. Other aspects of 
the behavior of 7r„ ooj such as zero distribution and interior asymptotics, were investigated 
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Figure 1. In this figure, K (respectively D) is the region enclosed by the 
black contour. Here, K has analytic boundary, and the dashed contour on 
the left is the inner-most contour outside of which we can find a univalent 
extension of 4>~^. Op is represented by the region outside of the contour 
corresponding to radius p. The curves outside of T are level lines of Pr- . 



in [15, 6, 7]. The following theorem provides an analog of [21, Theorem 1.2] for potential 
theoretic polynomials 7r„ g. 

Theorem 1. Let T — dK he a Jordan curve of class C'p+i."^ p + a> 1/2, and {t^u.s^^Zq^ 
be a sequence of polynomials satisfying orthogonality relations (9). Then, as n,s ^ oo, the 
leading coefficient Xn^s of TTn^s satisfies 



(13) 



ln+1 



Ik 



1 - 



n + l 



l + O 



1 



„2(p+a) 



Moreover, if T is an analytic Jordan curve, then the error terms in (13) can be replaced by 
0{p'^"-) for any p{T) < p < 1. It also holds that 



(14) 



'71+1 



1 - 



n + l 



$"$'[1 + (!]„)] 



uniformly on O as n, s ^ oo, where E„ is given by Table 1. 



Remark 1.1. When p — Q and limsup„ s^ao n/ s = \ the authors were unable to show that 
S„ — > as n — > oo (the employed method yields E„ = n^'^~"'), which is the reason this 
case is not included in Table 1. 

In general, the location of the zeros of 7r„^s depends on s (as well as, obviously, K and 
n). However, as the following proposition shows, this is not the case for a family of ellipses 
which interpolate between the unit circle and the interval [—2, 2]. 
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Analytic 




T 


p{T) <p<l 




lim sup n/s < 1 

n,s— >cxD 


limsupn/s = 1 

n,s— ^oo 




p= 1 


p = Q 


p = i 






log n 


log n 







Table 1. The error term I]„ depending on the smoothness of T. 



Proposition 2. Let q e [0, 1) and define (f){w) $ ^(ui) = u> + ^ so that K is the ellipse 
bounded by 0(T). Then, for all n < [s — 2J and all s including s — oo, 



ln+1 



n + 1 



/ I 



^2rl+2 



s — n — 1 
s + n + l 



1 



9 



<J)2n+2 



T/iat is, t/ie polynomials 7r„^s are i/ie renormalized Chebyshev polynomials of the second kind 
for the interval [—2y/q,2y/q'\, where ±2^^ are the foci of T . 

Remark 2.1. In the proposition above all the ellipses have unit logarithmic capacity (i.e., 
IK = !)• 



Remark 2.2. If g = 0, then T = T (/-iT = D) and 

/n+1 



Remark 2.3. If 9 = 1, then T = K = [-2, 2] and 

1 



1 - 



n + 1 



(15) -ninA^) 



s2 - (n + 1)2 



27rs 



Tl+l 



z - Vz^ - 4 



n+1 



It is easy to see that the asymptotic behavior of the normalizing constant in (15) is different 
from the one in (13). However, the case K ~ [—2, 2] is not covered by Theorem 1. 

Theorem 1 is the essential building block in proving results on asymptotic behavior of 
kernels K^^s and Kn,s defined in (10) and (11), respectively. 



Proposition 3. Let N < [s — IJ . Under the conditions of Theorem 1, it holds that 



Kn,s{z,w) 



N+1 



-{N + 1) 



$(z)<I>(w) 



N 



1 - 



$(z)$(w) 



Af+1 - 



1 - $(z)$(u;) 



1 - $(z)$(w;) 



(A^ + 2) 



N+1 



<I>(z)$(w) 



(16) 



1 - $(z)<I>(w) 
'O (max{l, AT^Sjv}^ 



1 - $(z)$(u>) 
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uniformly for z,w £ O , z ^ w, dist(z, dK) < const. /A^ and dist(w, dK) < const. /N. 

Moreover, it holds that 

(17) 

|$'(z)|2 riV(iV + l) / A^ + l\ N{N + 1){N + 2)' 



Kn,s{z,z) 

TT 

uniformly for z G dK . 



1 - 



6s 



-O (max{l, AT^Ea,}) 



To continue, denote by Aj^ the Hilbert space of holomorphic functions on D whose moduli 
are square-integrable with respect to the area measure. We equip with the norm induced 
by the inner product 



(18) {f,g) := / f{z)g{z)dA. 

J D 

Denote by K£){z,w), z,w E D, the reproducing kernel^ for Aj^. That is, 

(19) f{z) = / f{w)KDiz,w)dA 



JD 

for any / e It is known [8, Theorem 1.5.2] that Km,oo is the reproducing kernel for the 
set of polynomials of degree at most — 1 in the sense of (19) and that 

(20) \Kd{z,w) - Km,oo{z,w)\ ^ as TV -)■ oo 
locally uniformly for z,w G D. 

Theorem 4. Under the conditions of Theorem 1, (20) holds with ifjv.oo replaced by K^^s 
as N.s-^ OD, N <'ls~l\. 

Remark 4.1. The original proof of Theorem 4 as devised by the authors used the full strength 
of Theorem 1 and therefore excluded the case p = and i = 1. The authors are grateful to 
the anonymous referee who pointed out that only the asymptotics of Carleman polynomials 
TTn.oo on O is needed to prove the theorem and thus allowing all p > and £ G [0, 1]. 

To describe the asymptotic behavior of kernels (10) and (11) near the boundary dK, it 
is convenient to introduce the following notation. Set 

i/oM.^2?:il^ and H.ir):^, ''^^-\^^^\ 

and define i/^ to be the convex combination, 

3 — 3£ I 

(21) ^^(^)^^332^^°(^)+332^^^^^^' ^^(O'l)- 

Note that the value at the origin for each of these functions is determined by taking a limit; 
that is, -fff(O) — 1 for all I G [0, 1]. The following theorem is an analog of [13, Theorem 2.1]. 



Theorem b. Let N < [s-l\, z E dK, T{a, z) := a$'(z)$(z), £ := lim Ns^'^ G [0, 1] and 
for i> 0, set 

exp { - Re (r (a, z)) /£}, Re (r (a, z)) > 0, 



(22) u}{a,z) 



1, otherwise. 



^K£){z,w) = — ^^'^' ^ (^) _ where tp is any conformal map from D onto D [8, §1.5]. 
(i-^(z)^(™)] 
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Under the conditions of Theorem 1, assuming p > when £ — 1, it holds that 

23 lim ^ = Hi[T{a,z)+T{b,z)] . 

Moreover, if £ > 0, then 

(24) lim KnA^+W^^+w) ^u;{a,zMb,z)Hi(T{a,z) + Hb;z)) , 
and if £ = 0, then 

(25) ^ 

^nA^ +§^^+7f) Hp (rja, z) + r(6, z)) Re(r(a, ^)), Re(r(&, z)) < 0; 

N,s^oa KmA^i^) 1 otherwise. 

The convergence in (23)-(25) is uniform for a, b in compact subsets of C and z £ dK . 

Remark 5.1. The argument of r(a, z) is equal to the angle between a and $(2;)/$'(z), the 
outward normal to T at z. 

Remark 5.2. As is clear from (11), the function oj{a,z) is designed to describe the limit of 
P^^ (z + a/N) as N,s oo. This limit depends on whether or not the points z + a/N belong 
to O for N large enough. The case Re(r(a, z)) = corresponds to the situation when the 
sequence {z + a/N} approaches z e dK tangentially to the boundary. This does not cause 
a problem in (22) as this function is continuous with respect to a. However, when £ = 
formula (22) cannot be used as the limit is described by a discontinuous function of a and 
the convexity of the boundary dK at z starts to play a role. 

Remark 5.3. Observe that by putting s = oo (that is, £ = 0), formulae (16)-(23) specialize 
to the asymptotic formulae obtained in [13] for Carleman polynomials. Notice also that 
when s — N + 1 (£ = 1), the first summands in (16)-(23) disappear and only the second 
ones remain. For general £, formulae (16)-(23) turn out to be convex combination of these 
two extreme cases. 



3. Three Models of Potential Theoretic Ensembles 

Before proceeding to the proofs of our main results, we will present three models, a matrix 
model, an electrostatic model and a polynomial model, whose joint density of eigenvalues, 
particles and roots coincide with the potential theoretic ensembles we are considering. 

3.1. Entropic Normal Matrix Ensembles. The entropy of a self-map T on a metric 
space X is a measure of how the distance between nearby points is stretched under iteration 
of T. In the case where Z is an iV x iV complex matrix acting on C^, the entropy of Z is 
given by 

N 

MZ) = ^logmax{l,|A„|}, 

n=l 

where Ai, A2, . . . , Xn are the eigenvalues of Z [23]. We may use this to create a probability 
measure on normal N x N complex matrices, which we will denote by MnIC). 

There exists a canonical measure on A/'Ar(C) induced by the standard metric on C^^^ 
and we may define a probability density with respect to this measure by writing 

P^(Z) ^ ^e-'^^(^\ 

Am 
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where is a normahzation constant and s > iV is a real number necessary so that the 
probabihty measure is actually finite. 

This probability measure on normal matrices induces a symmetric probability measure 
on C''^ as identified with vectors of eigenvalues. This measure is absolutely continuous with 
respect to Lebesgue measure and its density is given as in (2) with w{\) — max{l, |A|}^^* 
[4, 16]. Normal matrix ensembles, and in particular the statistics of their eigenvalues, were 
first considered in [3] and [4]. 

The function A i-> logmaxjl, |A|} is the logarithmic (equilibrium) potential of the closed 
unit disk, and the weight for the entropic ensemble is formed from this in the obvious man- 
ner. We therefore see that the eigenvalue statistics of the entropic normal matrix ensemble 
coincides with the potential theoretic ensemble with K = D. 

3.2. Two-Dimensional Electrostatics. In two-dimensional electrostatics, charged parti- 
cles are identified with points in the extended complex plane. The potential energy of a 
system of two like charged particles located at z,w S C is proportional to — log|z — w\. 
More generally, if Zi, Z2, ■ ■ ■ , zn are the locations of A'' identically charged particles, then z 
determines the state of the system and the potential energy of this state is given by 

- X! 1^" ~ 

The energy is minimized when the particles are all at cx). In order for the system to be 
found in a state where the particles are at finite positions, there needs to be a potential (or 
other obstructions) which repels the particles from oo. We represent this field by V so that 
the interaction energy between a particle located at z and the field is given by V{z). The 
total potential energy of the system comprised of the N particles in the field is given by 

N 

E{z) = ^ V{zn) - ^ log|z„ - 

71—1 'm<n 

The system is assumed to be in contact with a heat reservoir so that the energy of the 
system is variable, but the temperature is fixed. In this setting, /3 denotes the reciprocal of 
the temperature, and the Boltzmann factor for the state z is given by 

^ n— 1 m<n 

This quantity gives the relative density of states, so that the probability (density) of finding 
the system in state z is given by 

^e-^^(^) where Zn ^ [ e-^^(=')dA^(z). 
Zn jc" 

Comparing with (2) we see that, when /? = 2 the density of states is identical with the 
density of eigenvalues of the normal matrix ensemble with weight 'w{z) = e~'^^^^K 

In this model, a compact set K is identified with a conducting region. A charge supported 
on K will distribute itself to minimize its potential energy, and this distribution, suitably 
normalized, leads to the equilibrium measure on K. In this way, we can think of the function 
—s\ogPK{z) as the potential energy felt by an oppositely charged particle at z when placed 
in the field given by the minimal energy configuration formed by placing a total charge of 
s on K. In this situation where our system consists of N charged particles, the condition 
that s > is required to make oo repulsive (or rather to make K sufficiently attractive so 
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that the particles do not flee to oo). It follows that the statistics of particles in this model 
agree with those of the potential theoretic ensemble for K. 

3.3. Roots of Random Polynomials. The Mahler measure of a polynomial f{x) £ C[x] 
is given by 

Af(/) = exp|^' log 1/(62-")! d^j, 

is an example of height function; that is a function which measures the complexity of arith- 
metic objects, in this case polynomials'^. One type of problem of interest to number theo- 
rists is to provide asymptotic estimates for the number of arithmetic objects whose height 
is bounded by C as C — )• c». For instance, for the Mahler measure, such estimates for the 
number of integer polynomials of fixed degree and Mahler measure bounded by C as C — cx) 
was given by Chern and Vaaler in [5]. They also gave a similar estimate for the number of 
polynomials with Gaussian integer (Z[i]) coefficients. 

In the latter case, the main term in their estimate came from the calculation of the 
Lebesgue measure of the set of polynomials of degree N with complex coefficients whose 
Mahler measure is at most 1. A key aspect of their proof is to show that this volume is 
equal to 

That is the volume is proportional to an integral of a (negative) power of the Mahler measure 
of monic polynomials with respect to Lebesgue measure on the non-leading coefficients of 
such polynomials. Moreover, after the change of variables from coefficients to roots of 
polynomials, this volume reduces to 
(26) 

IvFTT / ITTcxpI /'log|a„-e2..e|^J I TT |a„ - a„,|2 dA^(a); s = N + l. 
N + lJcN ]^^^^^ [Jo J J ^<^„ 

That is, this volume, up to the factor of tt/{N + 1) is equal to the normalization constant 
Zn for the potential theoretic ensemble for the unit circle for the value s = N + 1. In fact, 
Chern and Vaaler were able to show that this normalization constant as a function of s is a 
rational function in s with poles at positive integers < N. This striking result can be seen 
as a consequence of determinantal nature of the correlation functions. 

The derivation of (26) shows that the roots of a polynomial chosen randomly from the 
volume of complex polynomials of degree N and Mahler measure at most 1 obey the same 
statistics as those of the potential theoretic ensemble for the disk. This gives a polynomial 
model for these statistics. 

The computation of the normalization constant of this polynomial model for potentials 
for certain other compact regions (in particular the ellipses considered in Proposition 2) is 
given in [19], while a more general treatment for more general potentials is given in [20]. 
The special case where the family of ellipses degenerates to the interval [—2,2] on the real 
axis, and its application to the estimation of counting reciprocal polynomials with bounded 
Mahler measure is given in [18] and [19]. 



^Traditionally Mahler measure is used as a height of polynomials in Q[x], or more generally Q[x]. However, 
there is no obstruction in defining it for polynomials in C[x] 
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4. Proofs 

To prove Theorem 1, we use the method of normal moments in which we rely on the 
results in [21, Ch. I]. We must therefore discuss Faber polynomials before proceeding to the 
proof of Theorem 1. We start by stating several auxiliary facts that will be useful later. 

4.1. Auxiliary Facts. If g is holomorphic in O and vanishes at oo, then for each r G 
[l,oo) the restriction g{rT), r G T := {\w\ = 1}, can be interpreted as the trace on T of 
g{r/w), w ^ D, which is holomorphic and vanishes at the origin. The latter implies that 
JjT''g{rT)\dT\ = for all integers fc > 0. If, in addition, ft, is a positive function on [l,oo) 

and 



g{'w)h{\w\) 
theorem yields that 



(27) / g{w)h{\w\)dA 



is integrable with respect to dA for some j > 0, then the Fubini-Tonelli 



r^+^h{r)dr = 0. 



g{rT)\dT\ 



Furthermore, since dA{w) = \^'{z)\ dA{z), where w — <f>(z), it holds that 



(28) / ^^{z)<i>'(z)G(z)h{\<S>{z)\)dA = 



for any G holomorphic in O and vanishing at oo by (27) applied with g = {G o 4>)(l)' , where 
4> is the inverse of $ (granted I'l'-'^'G/id^l) | is integrable with respect to dA). 

In another connection, the Cauchy-Green identity for the domain D [12, Thm. 1.2.1] says 
that 

1 



(29) g{.z)h'{z) dA=-j^ g{z)Hz) dz 

whenever g and h' are holomorphic functions in D that continuously extend to T, where § 
always means integration in the counter-clockwise direction unless specified otherwise. Now, 
assume that g and h are holomorphic functions in O such that g has at least a double zero 
at infinity, and g, h and h' continuously extend to T. Then by using the transformation 
z H> 1/z and (29), one can show that the Cauchy-Green identity for O assumes the form 

1 



(30) 9{z)h'{z) dA=--j^ 9{z)Kz) dz. 

4.2. Faber Polynomials. Denote by Fn the n-th Faber polynomial for D associated with 
That is, 

J'j' t — z Ztti 

In other words, F„ is the polynomial part of Then it follows from Plemelj-Sokhotski 

formulae [9] that 

(31) F„ = $"$' + ^„ in O, 

where En is a holomorphic function vanishing at infinity with integral representation 

, . f $"(*)*'(*) dt 
Jrp t — z Zm 

We would like to point out that 
(33) En{z) ^o(\ 
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for all integers n > 0. Indeed, consider Fn+i, the {n + l)th Faber polynomial associated 
with 1. In this case (31) gets replaced by Fn+i = (f>"+^ +£'„+i, where En+i has an integral 
representation similar to (32). By differentiating both sides of the last equality, we get that 
(n + l)Fn = -F^^-i and (n+ !)£'„ = E'^^-^- As E^^i is holomorphic and vanishing at infinity. 
En has at least a double zero there. 

We are interested in the asymptotic behavior of 

(34) ml^:^ [ F,F\PK^'dA= [ FjFldA+ [ F,F\m-^'dA, j,fc<[s-2j, 

JC JD Jo 

where we used (12) for the second representation. 

It was shown in [21, Equation (1.32) combined with (30) above] that the first integral on 
the right-hand side of (34) can be written as 

(35) ^ - ^ I E,E-,dA) ^: ^ (4. + /.) , 

where the integral over O is well defined (finite) by (33). Moreover, it was also obtained 
there, see [21, Equation (1.45) and Lemma 1.5], that 

(36) \Id\<t-- ^"''.f' ^ or <const.p^■+^ 

where both constants are independent of j and fc, but depend on T and p (in the analytic 
case). Hereafter, by stating a double estimate of the form (36), we always assume that the 
first bound is given for T of class (7^+^'" and the second one for T analytic with p{T) < p < 1. 

On the other hand, the second integral on the right-hand side of (34) can be written with 
the help of (31) as 
(37) 



$i$fe|$'|2|$|"2s^^_^_ / $i$'£;^^|$|-2s^^_^ / £;^.$fc$'|$|-2^d^+ / EjEk\^\-'^''dA. 
Jo Jo Jo 

It can be immediately computed by conformality of $ that the first integral in (37) is equal 
to 

f $J$fe|$'|2|$|-2'*d^ = / w^^\w\-^'dA= ^ -4,-. 

Jo Jo s - (fc + 1) "-^ 

The second integral in (37) is zero by (28) applied with G — Ek and h{r) — r^^'', r E [1, oo). 
Conjugating the third integral in (37), analogous reasoning shows it is zero as well. Thus, 
similar to (35), the second integral on the right-hand side of (34) can be written as 

We claim that |/o| satisfies (36) as well, namely, 

(39) \Io\<^. ^ or |/o| <const.p^+'=. 

^ ' ' ' - (j + l)P+"(fc + 1)P+" I - 1^ 

Indeed, to prove (36), it was shown in [21, (1.44) and the following paragraph, (1.45)] that 

(40) ^|(£;„o0)(r)0'(r)|'|dr|<^-^^P^ or < const.p^" 

for all integers n > 0, where (/) is the inverse of $. Moreover, the monotonicity of L^- 
norms implies that if the integrand in (40) is evaluated at rr rather than at r for any fixed 
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r S (1, oo), the estimate remains valid. Then on accomit of 



< 



\{{E,Ek)ocl>){rT)\ \cb'{rT)\''\dT\ 



and the Cauchy-Schwarz inequahty, \Io\ is bounded by 



s- (fc + 1) 



|(i?fco0)(r)0'(r)hdr| 



1/2 



\{E,ocj,){T)cj,'{r)Y\dT\ 



1/2 



Clearly, (39) follows now from (40). 

Finally, gathering together (35) and (38), we get that 



s — STT (X s \ 



(41) 
where 

(42) ffe,, ID H W 



1 - 



k + l 



EjE, (1 - |<i>| 



dA 



const. 



and 

(43) kfc i\ < T- — TT—r, tvt;- 

^ ' (j + 1)P+" (A; + 1)P+" 

by (36) and (39). 



or |e^.j | < const. 



i+fc 



4.3. The Von Koch-Riesz Algebra. Denote by T) the algebra of all operators defined on 
lii^) by matrices A — \ak,j\j'k=^ respect to the standard basis for which 

1/2^ 

|A||i, :=max<| ^|afc^fe|, I ^'^W^j'^ 



< oo. 



k=0 



.j=Q fc=0 



It is known [11, Theorem II. 2.1] that if {A„} is a sequence in V converging to A e P (with 
respect to |j • ||p), then the determinant of I + A„ (I being the identity operator) converges 
to the determinant of I + A. 

Let {s„} be an increasing sequence of positive reals such that s„ — >■ oo as n — > oo. Set, for 
convenience, e^"^ :— when either j or k is greater than [s„ — 2J and define E,,^^ :— [el"j]ff^^Q- 
We also set Eqo ■— [^'kj]j°k=0' where we put 



(44) 



k + l 



EjEk dA. 



Observe that the estimate in (43) is also valid for s = oo. Using this bound, it is simple to 
verify that 

const. 



i|E,Jh < const. or ||E,J|p < - 



fe=i 



for each n including the case n = oo, where the constant YlV=i fc"^*^^"'""-' is finite as p + a > 
1/2. Thus, all the operators E^^ belong to the Von Koch-Riesz algebra V. Moreover, it 
holds that 



(45) 



IE., 



E, 



IX) 



as n oo. 
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Indeed, let {fc„} be a non-decreasing sequence of integers such that fc„ — >■ c» and fc„/n — >■ 

as n — oo. Then 

(46) 

const. 



- (fc„ + 1)?'+/^ " 



oo 



< const. p 



2k„ 



by (43). Furthermore, we can readily deduce from (42) and (44) using the notation of (35) 
and (38) that 



'k,k 



oo 



— \Io-Id\<'^ 



const. 



'k,k 



n n {k + l)2(p+") 

by (36) and (39) for all fc e {0, . . . , fc„ - 1}. Therefore, it holds that 

fe„-i 



fe.fc 



< — const, p 

n 



2k 



^k,k 



(47) E 

k=0 

Combining (46) and (47), we deduce that 

oo 

(48) E ' 



< const. 



^k,k ~ "=fc,fe 



-^0 as n — > oo 



by the choice of the sequence {fc„}. Analogously, one can show that 

2 

oo 



OO oo 



(49) 



as n — 7> oo. 



j:k = 



and therefore 



EE|<:^ 

j=0 k=0 

which finishes the proof of (45). 

Naturally [11, Section I.l], it holds that dct(I+E5^J = det Okj 
we deduce from the remark made at the beginning of this section that 
(50) det(I + Es,J ^ det(I + Eoo) > as n ^ oo, 

where the last inequality was shown in [21, Section 1.4]. 

4.4. Proof of Theorem 1. Since {Fn} is a complete system of polynomials, each 7r„^s can 
be expressed as a linear combination of i^oi ■ • • ? with the coefficients determined via the 
orthogonality relations (9). In fact, it holds that 



(51) 



7r„,s(z) = 



\/ Dn — l,sDn,s 



'n-l.O 



'1,1 



Foiz) F,{z) 



'"■n-l.n 
Fn{z) 



where the moments j are defined in (34) and -Dn,s ■— det 

Set A„ s ■— det[5kj + j]Jk=o observe that 
(52) 

n n 

t=o (^ + ^'^'^ - (fc + 1)) k=o k^] (^ + ''(•'* " + 1)) 

by (41), where the determinants Dn^sij) and A„ ^(j) are obtained from the same matrices 
as Dn,s and A„^s only with the last row and the (j + l)-st column removed. Given (43), it is 
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a straightforward algebraic computation using Hadamard's inequality, see [21, Lemma 1.7], 
to derive that 

const 

(53) An,.(j)< or A„,,0-) < consty+" 

for any j e {0, 1, . . . , n - 1}. 

On the other hand, the family {An^^} is bounded away from zero. Indeed, as mentioned 
just before (12), $'(oo) = "f'^^ and therefore the leading coefficient of is equal to 7^"~^. 
Hence, we get from (51) and (52) that 

/r-.N n+1 Dn-i.s n + 1 / n + l^A„„i.. 

(54) ,<„_,7^+i = ^— ^. — ^ 1 



Recall that any monic orthogonal polynomial has the smallest i^-norm with respect to the 
weight of orthogonality among all monic polynomials of the same degree. In particular. 



2n+2 s 

"'•n,n- 



Therefore, it follows from (54), (41), and (42) that 

A " + 1/^1 " + A„_i,, ^n + lf^ n+l\ , 

= (l + <.„)A„_i,, < A„_i,,. 
Hence, it holds that 

(55) inf min A„ s = inf Ai s_2l s > 

s l<ji<Ls-2j ' s L J. 

by (50) since A|^s_2j^s = det(I + E^), which proves the claim. 

Thus, expanding the determinant A„_s by the last row, we get that 

n-l 

(56) A„,, = (1 + e^, JA„_i,, + 5](-l)"+^-6,^,^.A„,,(j-). 

Dividing both sides of the equality above by A„_i.s and using (43), (53), and (55) yields 

Clearly, we get (13) by taking the reciprocal of (57) and substituting it into (54). 
Now, expanding the determinant in (51) by the last row as in (56) yields 

( p I -I y»+j- (j + 1)('S - J - 1) A„^,(j) \ 

[ h (n + l)(s-n-l) A„„i^/^J 



hi + 1 n+l\ I A, 



An 

Hence, by factoring out and using (31), the error term in (14) can be written as 



En ..„+,- (.7 + l)(l-^) A„,.(j) / ^\ ^ 

(n + 1)(1 - ^) A„_i,, V 
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Since |<i>| > 1 and |$'| is bounded away from zero in O, we get from (40), (53), (55), and 
(57) that the error term in (14) is of order 

iog(n+i) , 1 ij + m-'^) 1 



(-581 I ^ V 

^ ' (n + 1)P+" (n+ 1)P+" ^ 



(n + l)P+" ' (n+ 1)P+" ^ (n + 1)(1 - ^i±i) (j + 1)P+"' 

If hmsup„ n/s < 1, then the fractions (1 — '^)/(l — ^^) ^■''c uniformly bounded 
above and it easily follows from (58) that the error term in (14) is of order 

(59) 1°^^" + ^) ■ ' 



(n + l)P+" (n+ l)2(p+")-i' 

where the first summand is larger for all p > 1 and the second one is larger when p = 0. 
Clearly, the estimate for the error term in the case of analytic curve, can be derived in a 
similar fashion. On the other hand, if limsup„ g^ao = 1' then we use the estimate 

n+1 n+1 1 

1 — > 1 - 



3 n + 2 n + 2' 

which is valid since n < [s — 2J . In this case (59) gets replaced by 

logjn + 1) 1 

(n + 1)P+" (n+ l)2(p+")-2' 

where the first summand is larger for all p > 2 and and the second one is larger when p = 1 
(we exclude p = as in this case the above bound grows as n^^^^")). Analogous estimate 
shows that the error term in (14) is of order p"" when T is an analytic curve. This finishes 
the proof of Theorem 1. □ 

4.5. Proof of Proposition 2. Let Un be the monic Chebyshev polynomial of the second 
kind for the interval [—2y/q, "^y/q] . That is. 

It can be readily checked that the inverse of <I> is indeed (j>{w) ~ w + q/w, $ is the conformal 
map of the complement of [— 2y^, 2y^] onto {w : \w\ > with positive derivative at 
infinity, and the level lines of $ are ellipses with foci ±2y^. 

Let us show that polynomials J7„ are orthogonal on D with respect to area measure. It 
follows from the Cauchy-Green identity (29) that 

2i{k + l) [ Un{z)z^dA = i" ($"(z)-g"+V$"+2(z))$'(z)^dz 

(60) = / (r" - (7"+Vr"+2) (rq + l/T)*=+i(ir = 

Jt 

for k < n, where we used the identity t — 1/f on T and the last equality is a consequence 
of the facts /.j, T^r^dr — for all j > —n and T~''^^^T^dT = for all j < n. 
In another connection, it holds that 



Ur,{z)z'^\^{z)\-^'dA = / ($"(2)-9"+V3'"+'(-z))(^V^'W)l^(^)r'1*'(^)l''^^ 



(w" - g"+Vu;"+2),/)fe(u;)0'(w)|u;|-2^dA. 
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It is easy to check using the expressions (piw) = w + q/w and w = tt, r £ [1, oo) and r G T, 
that the chain of equahties above can be continued as 



(61) 



(rr)" - 



(rr) 



n+2 



1 - 



\dr\ 



rdr 



= 



' r qr^ 
r 

for k < n, since the Laurent polynomials in r integrated over T does not contain a constant 
term as the exponents of r range from n + k + 2 down to n — k and then from k — n down 
to —k — n. 

Altogether, the polynomials Un are orthogonal over C with respect to the measure 
P^'^^'dA. In fact, it can be easily shown that they are also the Faber polynomials for 
this K. It remains to compute the normalizing factor. Evaluating as in (60) and (61), we 
get that 



C/„(z)z"dA 



D 



1 1 

2i n + I Jf 



w 



,,n+2 



[qiu) 



ri+l 



1 



,,n+l 



dw 



and 



Un{zWMz)\~''dA 



(rr)" 



(rT)"+2 
„2n+2 



\dr\ 



(n+1) s + n + 1 



Thus, we deduce that 

-2 _ 



\Un{z)\'PK""iz)dA= I Un{z)z-^dA + / C/„ (z)z» | $(z) | "^^dA 



{n + l)(s - n - 1) 



D 



2n+2 • 



n — 1 



□ 



s + n + 1 ^ 

4.6. Proof of Proposition 3. Fix < [s — Ij and let C be a point such that 
(62) CeO and dist(C, K) < c/N 

for some fixed constant c. Further, let (q G T he such that |C ~ Col — dist(C, AT). Since 
|<I>(Co)| = 1 and is continuously differentiable in O (since T is at least C^'"-smooth) it 
holds that 

i<f (01 < 1 + mc) - *(Co)i < 1 + one - Col) - 1 + oiN-'), 

where the estimate 0{-) does not depend on the choice of C satisfying (62). Hence, 



(63) 



for some absolute constant. 



max 1$ (C)l < const. 
ke{i,...,N} 



Select z, w satisfying (62) and assume that z ^ w. Put, for brevity, u := ^{z)^{w). Then 
we get from the definition of ATat.s, (14), and (63) that 



Kn^s{z,w) = 



<I>'(2)$'(u;) 



$'(z)$'(u;) 



Af-l 



n=0 
/N-1 



N-l 



(64) 



n=0 
2s 



n=0 



-O (max{l,iV2l]jv}) ■ 
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Since, 

Y.{n + l)u- = -{N + l)- 



N-l M ^ _ ^N+1 



n=0 

and 



— u (1 — m)2 ' 



Y: in + Ifu- = -{N + lf-^ + iN+ 1) - + 2)7f%r + 2 VS^' 

1 — u (1 — uj^ (1 — U)^ (1 ~ U)'^ 



n=0 



the validity of (16) follows. In a similar but simpler fashion, we also get (17). 

4.7. Proof of Theorem 4. Recall the Christoffel variational principle: 

b(z)P 

(65) Kn^s{z,z)^ max ,^„-2s,a ^ z e C, 

dcg(p)<N \p\^Pj^^'dA 

and the reproducing property: 

(66) piz)= [ p{w)KNA^,^)PK^"HdA, deg(p)<iV, 

Jc 

that hold for all s G (1, oo] and < [s — IJ . It can be readily deduced from (65) that 

(67) Kn^s{z,z) < Kn,oo{z:,z) < Kd{z,z), 



where the second inequality follows from the fact that K]j{z,w) — X]^o ''^".oo('^)'''".oo('^) 
[8, § 1.5]. Furthermore, (66) together with (67) yield 

(68) / \KNA^^uj)\^dA< KDiz,z) and / \KN.siz,w)\\^{w)\-^''dA < Koiz, z). 

JD ' Jo 

It follows from (19) and (68) that 

\Kd{u,w)~ KN,s{'^^w)\'^dA = Kd{w,w)-2Kn.s{w,w)+ \KNju,w)\^dA 

D JD 

Therefore, (19), the Cauchy-Schwarz inequality, and the above estimate yield 

\Kd{z,w) - Kn,s{z,w)\ < / \Kd{u,w) - KN,s{u,w)\\KD{z,u)\dA 



JD 

< {KD{z,z)y^^{KD{w,w) - KN,si'w,w)y^^ . 

That is, we only need to demonstrate the convergence in (20) along the diagonal. Moreover, 
since (20) is valid for Kn,oo, it suffices to show only that Km,oo{z^ z) — Kj^^siz, z) — >■ locally 
uniformly in Z? as iV, s — > oo . To this end, observe that 

Kn,oo{z,z) = I KN,oo{z,u)KM,s{'^,z)Pj^'^%u)dA 

t^N.oo^^, i^J-t^N.syu, 

o 



KnAz,z)+ / KN.c<,{z,u)KNAu:z)\^{u)\-'^'dA 
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by (66) and the positivity of Kn^oo{z, z). Then we get from the equahty above, the Cauchy- 
Schwarz inequahty, and (68) that 

(69) < {Ku{z^z)) (^y^^^-^dAj . 

To esthnate the integral in (69), observe that 

|-f^jv,oo(2,u)P < ^JV,oo(^,2)-f^JV,oo(w, ") < Kd{z,z)Kn^^{u,u) 
by the Cauchy-Schwarz inequahty and (67). Hence, 

N-l 



I \KM,oo[z,u)\^mu)r^'dA < KD{z,z)y] f |7r„.o,n$|-2-d^ 
Jo ^^^Jo 

< const.KD{z,z)y] / |$|2"-2s|^/|2^^ 

n Jo 



n=Q 

N-1 



I \ n + 1 
< const.A£)(z, z) y 



s - (n + 1) 

n=0 ^ ' 

(70) = const./-i:i5(z,z) -, 

s ~ N 

where we used (14) for s = 00. Clearly, (70), (69), and the reduction process carried out 
above, prove Theorem 4 under the condition iV^/s — >■ as iV, s — > 00. The proof for all N 
readily follows from the obvious inequlity Km,s{z, z) < Km,s{z, z), N < M. □ 

4.8. Proof of Theorem 5. Since p > when f = 1, it holds that Sat — >■ as — >■ 00. 

Thus, we deduce from (17) that 

(71) i„ K„..i. ,in^fi_l,l\ 1^-^21 

N,s^oo TT \ 2 0/ TT 6 

uniformly for z ^ T. Fix c > and let a,b be such that |a|, \b\ < c. It follows from the 
Cauchy-Schwarz inequality that 



(72) 



a b 

Kn^,[z+-,z+- 



<max|ifAr,3(C,C)l, 



where the maximum is taken over all ^ satisfying (62). Furthermore, we get from the 
Bernstein- Walsh inequality and (63) that 

(73) \KnAC,0\ < |$(C)l'^'^"'Vax|i^w,,(w,w)| < const, max |AV.(«i,w)| 

for any ^ satisfying (62) with some absolute constant. Combining (72) and (73) with (71), we 
see that {K^^s {z + -f^^z + j^) A^~^} is a normal family for |a|, |6| < c, where the functions 
in this family are indexed by z e T, G N, and ,s G [N + l,oo). Therefore, it suffices to 
prove (23) only for those a, b for which z-|-^,z-|--^ ^ K (since T has a tangent at z, it holds 
that either z + a/N E K for all N large enough or z + a/A^ ^ iiT for all iV large enough). 
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As |$'| is bounded above in O, there exists a path 7 C O connecting z and ^ whose 
length is proportional to 1/A^. Hence, 



(74) 



a 
N 



$(z) - ^^'(z) = J ($'(0 - ^'iz))dt = « (^) , 



where the estimate holds uniformly for z G T and locally uniformly for a,b € 
|<i>(z)| = 1. This means that 



As 



(75) 
and 

(76) 



N 



b 

z H 

N 



1 



T{a, z) + t(6, z) 



o I — 

AT 



$ z 



a 



$ z + 



A^ 



N 



= exp 



{r(a,z)+T(6,z) + o(l)} 



where r(-, •) was defined before (22) and o(l) is again uniform for z € T and a, h in compact 
subsets of C. As before, we can assume without loss of generality that T(a, z) + t(6, z) 7^ 0. 
Then we get from (16), (75), (76), and the continuity of $' that 

6 



(77) lim Kn.s I z + 



|$'(z)p 



The limit in (23) now follows from 



uniformly for z g T and a, 5 in compact subsets of 
(71), (77), and (21). 

As obvious from (11), to prove (24), it suffices to show that 



(78) 



JV,, 



lim Pj^%z + a/N) = io{a,z) 



uniformly for z £ T and locally uniformly for a G C. To this end observe that an outward 
normal to T at z is given by $(z)/<I>'(z). Hence, the angle between the vectors a/N and 
$(z)/$'(z) is less than 7r/2 if and only if the vector a$'(z)/<i>(z) = T(a, z) belongs to the 
right half-plane. That is, if Re(T(a, z)) > 0. Hence, the hmit in (78) holds for Re(T(a, z)) < 
as Pk{z + a/N) = 1 for such a. Moreover, when Re(T(a, z)) > 0, (78) follows immediately 
from (76). The case Re(T(a, z)) = can be deduced by continuity and the uniformity of 
the estimate follows from the uniform character of the estimate in (74). Finally, the same 
arguments yield (24) for £ = 0. □ 



Appendix A. Plots of Correlation Functions 

To provide intuition for the results reported here we consider the scaled limit of Ri and 
i?2 of the entropic (potential theoretic, with K = 3) ensemble in a neighborhood of a point 
on the unit circle. By the radial symmetry of the weight it suffices to restrict ourselves to 
a neighborhood of 1. In this case T(a, 1) = a and uj{a,l) = g-^'^i'^')/^ for Re(a) > and 
uj{a, 1) = 1 otherwise. As before, if £ = and Re(a) > then we take a;(a, 1) — 0. For 
convenience we define the scaled kernel at z = 1 by 

Hi{a, b) = uj{a, l)uj{b, \)Hi (a + 6) . 

The limiting density of scaled eigenvalues is then given by 

R[{a) = Hi{a,a), 
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and the scaling limit of the second correlation function is given by 

Ri{a, b) = Hiia, a)He{b, h) - Hi{a, b)He{b, a) 

The visualizations provided here are for the cases where a and b are either real or on the 
imaginary axis. 

A.l. Tangent to the curve. The tangent line of the circle at z = 1 is parallel to the 
imaginary axis and the local density of eigenvalues in this direction is given by 

Ht{it,it) = 1. 

This is expected since the spatial density of eigenvalues on the unit circle must be invariant 
under rotation (and locally, this rotation is given by translation up the imaginary axis). 

Looking at the second correlation function, when a and b are on the imaginary axis, we 
see that is a function oi t — —i{a + b). Figures 2, 3 and 4 show plots of the second 
correlation function for various values of £ in various regions as a function of t. 




Figure 2. plots of R\{a, b) as a function oit = —i{a + b). The second plot 
is an enlargement of the shaded region. 




Figure 3. plots of i?§(a, b) as a function of t = —i{a + b). 

By way of comparison we also provide plots of the second scaled correlation function for 
ensembles with the sine kernel. Specifically, 

sin iia-b) 12) 

S(a,b) = 2 V rr^ and RT(a,b) = 1 - S(a,bf . 

[a-b) 
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(The slightly unusual normalization given by the superfluous appearing factors of 2 in the 
first equation arises in the scaling limit when we take the expected distance between eigen- 
values to be 27r — this allows for the most accurate comparison with our other figures). 




5 10 15 20 25 30 ' 5 10 15 20 25 30 



Figure 5. plots of -R2"(a, b) as a function of i = a — 6. 



A. 2. Normal to the curve. In the regime where a and b are real, we are looking in a 
neighborhood of z = 1 in a direction perpendicular to that where the density of eigenvalues 
becomes constant. That is, the first scaled correlation function should decay as a moves 
away from 0. Negative a corresponds to moving into K (where the potential is constant) 
whereas positive a corresponds to moving away from K where the potential acts to make 
oo repulsive. As i decreases to 0, the field increases in strength until at £ = there is no 
possibility that an eigenvalue can be outside K. That is, when £ = the first correlation 
function vanishes for a > 0. 
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Figure 6. Plot of the interpolation between R\{a) and R\{a) as a function 
oit — — Re(a). When £ = there is a sharp cutoff at t = 0. 

When a and h are real, R2{a,b) is no longer a function of a linear combination of a and 
6, and we plot this as a surface for £ = and £ = 1. 




a 



Figure 7. Plot of R\{a, h) when a and b are real, with part of the surface 
removed to see the cross-section. 
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